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ABSTRACT, There are many ways of programming an actuator or effector 
for movement between th same two points. In the interest of efficiency it 
is sometimes desirable to program that trajectory which requires the least 
amount of energy. This paper considers the minimum energy movement for 
a spring-like actuator abstracted from muscle mechanics and energetics. It 
is proved that for this actuator a bang-coast-bang actuation pattern 
minimizes the energy expenditure. For some parameter values this pattern 
is modified by a singular arc at the first switching point. A surprising 
limitation on the duration of coast is demonstrated. Some relaxations of the 
restrictions underlying the spring model are shown to preserve the bang- 
coast-bang solution. 
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Though considerable effort has been expended in the study of the 
human motor system, the execution of even simple movements is not Mel I 
understood. One current theory holds that movements are memorized in terms 
of final position [Bizzi et alJ. The organism selects length-tension 
curves of agonist and antagonist by adjusting the innervation levels so 
that the intersection of those curves occurs at the desired position 
[Feldman], 

This process is illustrated by the hypothetical length-tension 
curves of agonist and antagonist muscles in figure 1. Suppose the system 
is currently at length L under innervation rates gj for the agonist and n| 
for the antagonist. If the innervation rate of the agonist is changed to 
g 3 , a different agonist length-tension curve is selected and the 
/*-*% equilibrium length shifts to Lj. Assuming no delay in tension development 

and ignoring velocity effects, the arrow in the figure indicates the 
tension course. There is an isometric buildup of tension from P to P 2 
followed by an isotonic decay to P Jt where the tension in agonist balances 
the tension in antagonist. The final position theory maintains that the 
position Lj can be reached independent of starting position merely by 
setting rates n t and g 3 . This theory is interesting from a manipulation 
viewpoint because it obviates the need for precise trajectory calculation. 

There are many choices of agonist-antagonist length-tension curve 
pairs that have Lj as equilibrium position. One choice that could be 
expected to require less energy is n 2 and g 2 , which minimize the isometric 
tensions. More generally, it is conceivable that some complex sequence of 
innervation rates {n it gj} might require less energy than a scheme which 




^0 "^i length 
Figure 1. The equilibrium point of the intersecting 
length-tension curves of agonist (g labels) and 
antagonist (n labels)' shifts from L to L when the 

firing rate of the agonist is raised from g to q 

^1 ^3 

and the antagonist rate remains at n . 
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Figure 2 . Length-tension curves from the cat soleus 
muscle (Rack and Westbury) . 
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select 8 the final length-tension curves immediately. The determinaton of 
this optimal innervation pattern is the focus of this paper. 

I. Muscle Mechanics 

The investigated properties of muscle present a too complicated 
view for analytic treatment. The plan in my research has been to simplify 
the muscle mechanics until an analytic solution to the optimal energy 
problem is possible, then to examine if the nature of the solution is 
changed by adding some of the excluded muscle properties. A full-blown 
formulation of the problem would require numerical methods for solution, 
and I intend eventually to carry out this analysis. 

The length-tension curves in figure 1 are themselves a 
simplification of actual length-tension curves such as in figure 2. The 
simplification arises from extracting only the linear portions at short 
lengths. There is some question as to what portion of the length-tension 
curves are used in actual movement. Some authors [Zierler, Hill 1970, Cook 
and Stark] maintain that the active portions occur near L , the length at 
which there is maximum isometric tension and which Zierler also calls the 
rest or natural length of muscle. Collins et al. [19751 on the other hand 
report that the linear portions are used in eye movement. 

For the present we accept the hypothetical length-tension curves of 
figure 1. The curves are also assumed parallel, as reported by Collins et 
al.; Rack and Uestbury, however, report a decrease in slope with firing 
rate for this linear portion (figure 2). The curves in figure 1 lead to a 
model of muscle as a spring with variable zero setting. The slope K of the 
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curves represents the spring constant, and the variable zero setting L z 
corresponds to the selection of firing rate* The force exerted by a muscle 
is thus K(LrL E ). 

An important simplification of muscle properties is to exclude the 
series elastic component. Ue also exclude the parallel elastic component 
and the active damping. The resulting muscle model is presented in figure 
3# The equation of motion for the spring system of figure 3 is: 

m dfx » - bdx + k f (x r x) - k n (x-x ) (1) 

dt 2 « 

Define a control variable X and a state variable U as be low. 
. t xj ] C x 1 [ u, ] t x. ] 

x-r ].[] u«t ]-[] 

t x 2 ] [ x ] [u 2 l (kJ 

Setting the mass m - l f the state variable representation of the spring 
system is: 

. [81] [80] 

X - t J X + I J U (2) 



[ -k r k n -b I C k s k n ] 



(lore compact I y t 



X - A X + B U (3) 

• f(X,M) 

II. Muscle Energetics 

The energy E expended during movement equals work plus heat. The 
work U may be subdivided into conservative work performed on the mass m and 
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Figure 3. Simplified muscle model with equation of motion: 



where 



m x = - b x + k (x -x) - k (x-x ) 

9 g n n 



b * coefficient of passive damping 

k g = spring constant of the agonist 

k n = spring constant of the antagonist 

m = mass 

x * position 

x g = variable zero setting for the agonist 

X n = variable zero setting for the antagonist 
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nonconservative work performed on the viscous element. The isometric heat 

Qj is given off in maintaining the muscle at a particular tension P . The 

rate of energy expenditure is thus 

• . 

E - Po v (power) + a P (maintenance heat rate) 

- (x 2 + a) (k g (u r xj) + k n (x,-u 2 )) (4) 

where v is velocity and cc is the maintenance heat coefficient. The two 
force terms have been summed because each contributes to energy loss. Ue 
have excluded the shortening heat because the active damping was also 
excluded, and because there may be a theoretical relationship between the 
two [Huxley, Cap I an]. The transient characteristics of heat production 
have also been excluded. 

III. The Euler-Lagrange Equations 

The task now is to find the time varying control U(t) that 
minimizes the energy used in moving between two points in a fixed interval 
of time. Let V represent the energy consumed in applying the control U to 
yield the trajectory X. The problem of minimizing V is readily approached 
by techniques of modern control theory. The fundamental equations that the 
optimal control U(t) must satisfy are derived from a theorem from the 
calculus of variations . This theorem states that in o^der to find the n- 
vector X(t) that minimizes V(X), where 



L(X,X,t) 



V(X) - / L(X,X,t) dt (5) 

subject to the constraint relations 

♦ 
gj(X,X,t)-8 i - l,...,m < n (6) 
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then X(t) satisfies the Euler equations 

BV (X.X.t) - d 3LMX.X.t) « 8 i«l n (7) 

3k j dt 3xj 

where 

LMX f X,t) -L(X,X,t) +S Xj(t)gj(X f X,t) (8) 

and Xj(t) f s are the multiplier functions tSchultz and Melsal. 

Applying this theorem to the optimal control problem, the state 
equations X ■ f(X,U f t) represent the equality constraints* L represents 
the rate of change of energy E. The Hamilton! an H « L + X^f represents 
(8), where X T - l\\ \zl • By applying the Euler equation first for X and 
then for U, it can be shown that the minimizing U(t) satisfies the 
following two Euler-Lagrange equations tSchultz and Me I sal. 

X - -H x (9) 

H u « 8 (18) 

IV. The Minimum Principle 

Because L is linear in the control U f there will not generally 
exist a minimum energy solution. To obtain a realistic solution, 
constraints must be placed on the control. The solution in this case will 
Me on the constraint boundaries EBryson and Hoi. Constraints on U f 
however, make it impossible to differentiate H with respect to U. 

The minimum principle of Pontryagin makes It possible to proceed 
from this point. Pontryagin showed that even if the control is 
constrained, one still obtains a minimal solution by finding the u° » 
u°(X f Xft) to minimize the Hamiltonian H f but by inspection rather than by 
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differentiation* After finding the minimizing u° f one forma H° » 

H(X t u° f X.t) and then solves the following tuo equations tSchuitz and 

lie I sal. 

X - 3j£ (11) 

3X 

X « - 3H° (12) 

3X 

There are two natural constraints that fall on the control U. 

First, the spring cannot push. 

U| - xj > (13) 

x, - u 2 > 8 (14) 

Second, springs have a maximum tension that they can exert. Without this 

constraint the solution would involve an infinite impulse. For the moment 

we assume the maximum tension is constant and independent of length: 

U| - X! < C| (15) 

Xj - u 2 £ c 2 (16) 

where cj and c 2 are constants. The case of maximum tension varying with 

length is deferred until section XI. 

V* A Bang-Coast-Bang Solution 

To facilitate inspection of the Hamiltonian, we expand H « X T f + L 
into three lines, the first depending on u lf the second on u 2 , and the 
third on neither control. 

H - k f (uj-Xi) («+x 2 +X 2 ) (17) 

+ k n (x,-u 2 )(a+x 2 -X 2 ) 
+ x 2 (X r bX 2 ) 
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To minimize H with respect to u lf we observe that if a+x 2 +X2 > 9 th * n H is 
minimized when u^xj. If a+x 2 +X2 < then H is minimized with Ui-xj+cj. 
Similarly, it can be shown for u 2 that when a+x 2 -X 2 < the minimizing u 2 
lies at xj-c 2 ; otherwise u 2 is at x t . Combining these results, one finds a 
bang-coast-bang solution to the minimum energy for muscle movement. 

Case 1; X 2 < -(ot+x 2 ) 

Then Uj«Xj+c |f u 2 «Xj. 
Case 2s |X 2 | £ («+x 2 ) 

Then Uj»X|, u 2 «Xj. 
Case 3t X2 > («+x 2 ) 

Then U|»x |f u 2 «x r c 2 . 

VI. The Solution Equations 

Substituting the minimizing u° into H, one obtains three functions 
corresponding to the three cases. 

Case Is VP « k g cj(a+x 2 +X 2 ) + x 2 (Xi-bX2) 

Case 2s H° - x 2 (Xi-bX 2 > 

Case 3: M° - k n c 2 (a+x 2 -X 2 ) + x 2 (X r bX 2 ) 

The differential equation (11) and its solution becomes for the three 
cases: 



, [01] [ ] 

X - [ 1 X+ I ] 

(0-b) I kc 1 
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x 2 (t) - x 2 (t )e~ b(t - to) + kc(l - e" fa{t ^ to) ) (18) 

b 

x,(t) - X|(t ) + l(x 2 (t ) - x 2 (t)) + kc (t-t ) (19) 

b b 

Mhere c«C| and k«k g for case 1; c-0 for case 2; and c«-c 2 and k«k n for 

case 3* The time t represents the starting time. The differential 

equation (12) and its solution are: 



[00] t 1 

X - I 1 X - t 1 

t -1 b 1 t kc J 



Xi(t) "- Xi"Ct ) (20) 

X 2 (t) - X2(to)e b(t "" to) + Xi(t ft )+kc (1 - e b(t " to) ) (21) 

b 

where c and k have the same meaning as above except c«+c 2 for case 3. 

Since X| Ct) is constant, it appears henceforth as X) without a time 

dependence. 

VII. The Extremal Versus Singular Solution 

It is proven in appendix A that there are exactly three events In 
the extremal bang-coast-bang solutions an acceleration period, a coast 
period, and a deceleration period. No other combination of bangs and 
coasts is minimizing. However, a nonextremal minimizing solution may arise 
from a singular arc at the switching points. The Hamiltonian (17) has the 
curious property that if X2-|«+x 2 | then the corresponding control may take 
on any value and still minimize H. If a control can be found to maintain 
X2— I «+x 2 1 for a finite time interval, then a non-extremal solution to 
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energy minimization might exist. This situation is called a singular arc 
and arises from a performance index linear in control but quadratic in 
state IBryson and Hoi. To maintain X2*l«+X2l for a finite time interval, 
all time derivatives of the two switching curves must be zero: 

d n tXrH>+x,) - 8 n*8 (221 

dt n 

d n (X»-«-x,) - 8 n*8 (23) 

dt" 

Carrying through the analysis for a singular arc at the first switching 
point (22) f the time varying force during the singular arc is* 

k g (u,-x,) • X, + ba + 3bx 2 (ti)e 2b(t ~ t|) ♦ 3(b«+Xi) (e 2b(t " tl> -l) (24) 

2 

Unfortunately no sufficient condition has yet been developed to 
test whether a singular arc is minimi zing 9 and one must compare values of 
the performance index for specific parameter values for the singular arc 
solution versus the extremal solution. Depending on the choice of Xj the 
force (24) takes one of the three forms in figure 4, 

Of these forms only 4C has been found minimizing for some parameter 
combinations. To search for such combinations, a set of parameters was 
initially deduced from Rack and Uestbury (table I). The elapsed distance 
x f and the elapsed time t f are variable and have been chosen as 8.2 cm and 
8.4 sec respectively. The initial and final velocities are assumed zero. 
For the extremal solution (bang, coast, bang) there result 8 nonlinear 
equations in 8 unknowns from (18) -(21) and the initial conditions fx f ,t f ). 
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time 



t2 



Figure^. The three possible forms for force in a singular 
arc solution. 
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Table I 

k « 2 kg/cm 

m - 0.2 kg (plausible value) 

b/m » 3.16/sec (chosen to give f«0.5) 
c ■ 1 cm 
a » 0.1 cm/sec (deduced from Uo I edge) 

For the singular solution (bang, 4C, coast, bang) IB nonlinear equations in 15 
unknowns result from (18)-(24) and the initial conditions. The equations 
were solved numerically by Newton-Raphson and gradient methods. Individual 
parameters were varied and energies of movement computed from (5). Solving 
(5), the energy for the extremal solution is 

E - kc(dx|(t|)+dx|(t 3 )) + kcaOt|+3t 3 ) (25) 

where t t is the switching time from acceleration to coast, t 3 is the time 
at the end of deceleration, 3t| is the duration of acceleration, 3x|(tj), Is 
the distance moved during acceleration, and dx t (t 3 ) is the distance moved 
during deceleration. For the singular solution, the energy is 



4) -f / k g (uj-Xj) (x 2 (t)- 
/ft 



E - kcOxjttjJ+ax^t*)) + kcccOt,+3t 4 ) + / k f (u r x,) (x 2 (t)+cc) dt (2B) 

where t| is the switching time from acceleration to the singular arc 4C f t 2 
is the switching time from 4C to coast, and t$ is the time at the end of 
deceleration. The force k g (u r xj) is given by (24), while the velocity 
x 2 (t) is 

m #* 1 2b(t-t,) u . , 2b(t-t,) ._ ,__ 

x 2 (t) m x 2 (tj)e l + bq+Xi (e ■ -1) (27) 

2b 
The energies for the extremal versus the singular solution are 
compared in tables Ila-g; the units are kg cm/kg wt. In table lie 
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Table Ha 




Table lib 




k Singular 


Extremal 


b 


Singular 


Extremal 


7 impossible 


0.708 


2.3 


impossible 


0.545 


8 0.B6U7 


0.66106 


2.4 


0.552036 


0.552034 


*10 0.6173 


0.6164 


2.6 


0.5675 


0.5674 


20 0.5640 


0.5627 


3.0 


0.6020 


0.6014 


30 0.5513 


0.5509 


*ie 


0.6173 


0.6164 


40 0.54554 


0.54552 


4.0 


0.7045 


0.7030 


50 0.5423 


0.5426 


5.0 


0.820 


0.824 


100 0.536 


0.537 


10.0 


1.45 


1.51 






15.0 


2.13248 


2.13246 






15.557 


impossible 


2.4 



Table He (b*ifk) 



Table lid 
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k 


Singular 


Extremal 


c 


Singular 


Extremal 


7 


impossible 


0.664 


0.7 


impossible 


0.708 


8 


0.632026 


0.632023 


0.8 


0.6612 


0.6611 


*10 


0.6173 


0.6164 


*1.0 


0.6173 


0.6164 


14 


0.6442 


0.6440 


2.0 


0.5640 


0.5627 


20 


0.708 


0.725 


3.0 


0.5513 


0.5508 








4.0 


0.54554 


0.54551 








5.0 


0.5427 


0.5437 








6.0 


0.5401 


0.5407 



Tabic 


» He 




Table IIf(ksi6. b=4) 




tf 


Singular 


Extremal 


x, 


Singular 


Extremal 


0.35 


impossible 


0.814 


0.006 


impossible 


0.0041 


0.36 


0.76123 


0.76122 


0.05 


0.0632 


0.0634 


*0.4 


0.6173 


0.6164 


0.1 


0.1925 


0.1942 


0.45 


0.512 


0.511 


*0.2 


0.6647 


0.6681 


0.5 


0.444 


0.445 


0.3 


1.4636 


1.4620 


0.6 


0.358 


0.369 


0.35 


2.0120 


2.0072 








0.4 


2.6878 


2.6825 








0.54 


impossible 


6.4598 
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Tabic IIg (k=16. bM) 



a 


Singular 


Extremal 


0.06 


8.6112 


0.6161 


♦8.1 


8.6647 


0.6681 


8.2 


8.7974 


0.7981 


8.25 


8.8633 


0.8631 


8.5 


1.1903 


1.1882 


1.8 


1.8393 


1.8382 


1.7 


2.748417 


2.748411 


1.8 


impossible 


2.8784 



the parameters k and b are varied simultaneously but at a fixed damping 
ratio of 0,5. In tables Ilf-g the parameters k and b are respectively set 
at 16 and 4 rather than at the table I values Mhere the extremal solution 
is minimizing over the whole range of x f and a. The initial values from 
table I are starred in table II. A singular solution becomes minimizing 
with high values of k, b, c, and t f , and with low values of a and x f . As 
the parameters cause the coast time to approach zero (higher b and x f , 
lower k f c f and t f ), the singular and extremal solutions become identical 
because the 4C portion vanishes. 

For the extremal solution it is proved in appendix B that there is 
an upper limit on the duration of coast. It is tempting to speculate that 
for longer coast durations a singular solution becomes minimizing, but the 
singular solution in table II is not always minimizing under these 
conditions. Perhaps a different combination of bangs f coasts, and singular 
arcs would then be minimizing, but this remains an open question. Some 
combinations can be proved impossible, such as (bang, 4C, coast, 4A-C, 
bang! . 
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VIII. Spring Model Relaxations: X Dependencies 

A natural question is whether the minimum energy solution is 
changed by incorporating a more realistic muscle model. For those 
relaxations of the spring model involving only X dependencies/ the answer 
is that the solution remains bang-coast-bang. The reason is that the 
Hamiltonian H remains linear in the control U, and the minimization of H 
with respect to U occurs at fixed X. Whether the solution also remains 
acceleration-coast-deceleration needs to be determined for each case. 

Relaxations of the spring model involving X dependencies include 
the following. 

1. Position Limits on Tension 

For real muscle the maximum isometric tension varies with position 
(figure 2). This makes cj and c 2 into functions of x lf but the controls 
will still fall at the extremes wherever they are 9 

2. Velocity Limits on Tension 

Actual muscle exhibits a hyperbolic force-velocity relation. If P 
is the isometric tension, then the maximum force P that can be produced for 
a velocity v is [Hill 19381 (see figure 5): 

P - Po - v tPn+a) (28) 

v+b* 

The term (P +a)/(v+b f ) can be considered the coefficient of active damping. 
The coefficient a has been determined as .2B P ; the force P then becomes 
P -1.25P v/(v+b , J. The literature conflicts on the value of active damping 
during lengthening. For consistency with the shortening heat (below), it 




Tension 



Figure 5, Tension dependence on velocity (Hill 1938) 
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Figure 6. Hypothetical length-tension curves with 
the property that at any given length the slopes are 

the same for all choices of u . 

2 
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is assumed the same as active damping during shortening. 

Associated Hith the active damping is an extra heat expenditure 
above the isometric heat due to shortening. This shortening heat rate is 
[Hill 19641: 

Q 8 - .IB P v + .18 Pv (29) 

The isometric heat rate Qj remains aP , but the power is now Pv. 
Substituting the sum of spring forces for the isometric tension P and (28) 
for P f the energy rate is* 

L - (k g (u,-x,) + k n (u 2 -xj))(a - B.09x 2 + 1.45b' ) (30) 

xg+b* 

Similarly it can be show that the equation of motion is: 

x 2 • - b x 2 + (1 - 1.25x,) (k t (ui-xi) - k n (x,-u 2 )) (31) 

X2+b f 

When these terms are combined to form the Hamiltonian. the control is seen 
to remain linear. Hence the solution is once again bang-coast-bang. 

3. Spring Constant Variations with Position 

One way of bringing the simplified length-tension curves of figure 
1 closer to those of figure 2 is i I lustrated in figure S. The spring 
constant k n varies with position, but at any given position the constant k n 
is the same for all controls u 2 . Under these conditions the solution 
remains bang-coast-bang. 

4. Parallel and Series Elastic Elements 

The incorporation of these elements into the model is depicted in 
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figure 7. Since the parallel elastic element depends only on position, it 
does not change the solution. The series elastic elements and the active 
springs may be replaced with equivalent springs with constants 
k f * *k f k 8 /(k f +k s ) and k n *«k n k 8 /(k n +k s ). This modification also has no effect 
on the solution. 

IX. Spring Model Relaxations: U Dependencies 

In figure 2 the spring constant k n is seen to vary with firing rate 
at any fixed position. The linear portions of these length-tension curves 
when extended seem to intersect at a common point (figure 8). In this 
circumstance the spring constant k„ is T /(u 2 -x ). The u 2 terms of the 
Hamilton! an H become: 

H f - T ft (xi-u») («+x 2 -X 2 ) 
u 2 -x 

- T («+x 2 -A 2 ) ( xi-x fl - l) (32) 

u 2 -x 

If «+x 2 -X 2 > then H f is minimized at u 2 «xj; otherwise u 2 «x r c 2 . That is 
to say, the solution for u 2 is exactly the same as in section V. A similar 
analysts holds for up Thus the minimizing pattern is also bang-coast- 
bang. 
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Figure 7. An expanded muscle model incorporating series 
elastic elements k and a parallel elastic element ~k . 




length 



Figure 8 . The length-tension curves of Rack and Westbury 

when extended meet at a point. The dependence of k on u 

n 2 
can be characterized by T /(u -x ) . 
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Appendix A 

In this appendix it is shown there are only two bangs and one coast 
in the extremal solutions one acceleration, followed by one coast period, 
terminated by one deceleration. No other combinations of coasts and bangs 
are possible. To demonstrate this is the only possible combination, it is 
necessary to examine the switching curves and their time derivatives. 

The first lemma shows that once the control has passed from 
acceleration to coast, then the control cannot return to another 
acceleration but must proceed to deceleration. The second lemma shows that 
once deceleration has started, the deceleration must continue until the end 
of the movement. This proves that the acceleration-coast-deceleration 
combination is the only possible one. In the following it is presumed that 
m*^ the movement starts with acceleration in the positive xj direction. Hence 

all velocities are positive. 

Lemma It After acceleration, the glide period cannot double back to 
another acceleration. 

Proof; The proof of this lemma proceeds by examining the time derivative 
of the acceleration-coast switching curve (henceforth referred to as the 
elope of the switching curve). The slope of this curve is initially 
positive at the transition from acceleration to coast. In order for 
another acceleration to follow the coast period, this slope must become 
negative, leading to a contradiction. 

At the first switching time t| the acceleration-coast switching 
curve is zero. 
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X 2 (t,) + a + x 2 (t,) - 8 (AD 

After the acceleration period, the coast equations are: 

% *±\ % # A % b(t-ti) #4 b(t-ti) % # A «* 

X2(t) » X 2 ttj)e l + 3^(1 - e ' ) (A2) 

b 



x 2 (t) - - (X 2 Ct l )4«)e" btt " t|) 



(A3) 



Thue 



X 2 (t)+«+x 2 (t) • X 2 (t,)e l + jLLll-a ' > (A4) 

b 

f% r. i , -b(t-tj) 
+ * - (X 2 (t|)+a)e ■ 

The slope of this switching curve is: 

e b(t ~ t,} (bX 2 (t,)-Xi) + b(X 2 (t,)+a)e" b - t " t,} (A5) 

At t-tj, the slope of the switching curve is: 

2bX 2 (tj)-Xi+ba > (A8) 

One can show this quantity cannot be less than zero. Next, suppose the 
coast doubles back to another acceleration. At some point the slope must 
go through zero. This time t is found from (A5) as: 

e 2b(t ~ t,) - b(X,(t t )+q) > 1 (A7) 

X r bX 2 (t,) 

Case 1: Xi-bX 2 (t|) > 0. 

Then X 2 (ti)+a>0, contradicting (Al). 
Case 2: Xj-bX 2 (t|) < 0. 

Crossmultiplying (A7) and col lecting terms, 

2bX 2 (t|)-Xi+ba<0, contradicting (A6). 
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Thus after acceleration, the coast period must eventually arrive at the 
deceleration switching point. 

Lemma 2: The movement is locked in deceleration until the end. 
Proof: It Hill be shown that if deceleration ever switches to coast, then 
the 8 1 ope of the coast-deceleration switching curve requires an immediate 
return to deceleration. Hence the movement is locked in deceleration until 
the end. 

Suppose there is a time t 3 when deceleration switches to coast. At 
this point the coast-deceleration switching curve is zero. 

X 2 <t 3 )-«-x 2 (t 3 ) - 8 (A8) 

The coast switching curve XgftJ-or-^tt) is 

* n x b(t-t 3 ) * M b(t-t 3 K , A , -b(t-t 3 ) 

X2<t 3 )e * + &j.(l-e ) - « - x 2 (t 3 )e 3 (A9) 

b 
The 8 1 ope of (A9) is: 

e b(t ~ t3) (bX 2 (t 3 )-A,) + bx 2 (t 3 )e~ b(t ~ t3) (A18) 

At time t 3 the slope (A18) is bX2(t 3 )-X)+bx 2 (t 3 ). This is positive since 
X2<t 3 )>8 f x 2 (t 3 )>8, and X|<8 (lemma 3). This means that deceleration would 
bounce off the coast boundary and immediately continue the deceleration. 
Furthermore* since the slope is positive* the deceleration would not 
immediately switch back to coasting, causing chattering. 
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Lemma 3t X| < 0. 

Proof: At the second switching point t 2 we have 

X 2 (t 2 )-o-x 2 (t 2 ) .0 (All) 

From (A2) and (A3), this becomes 

* ii i b(t 2 -ti) . ,, b(t 2 -ti) -b(t 2 -ti) ,»,»»» 

X 2 (t|)e c ' + ij.d-6 * ' ) - « + (X 2 (t,)+«)e z ' . (A12) 

b 
Rearranging, 

. .. b(t 2 -t,), . (lll b(Vt|) -b(t 2 -t,). ., -b(t 2 -t,). ..... 

^(1-6 c ' ) - - \ 2 (t,)(e *■ ' +e z ' ) + «(l-e z ' ) (A13) 

b 
From (Al) and (A3) we find an expression for X 2 (t). 

A 2 (t,) - - a - k l c 1 (l-e" b(trt ° ) ) (A14) 

Substituting into (A13), 

% ii b(t 2-*lK ,. „ -b(t,-t ) w b(t 2 -t,) -b(t 2 -t,) ■ ,_,__ 
fcj.(l-e * * ) - k.ci(l-e ' ° ) (e z ' +e z ' ) ) (A15) 

b b 

+ «(l+e ' ' ) 



Thus 

x, . K.c,a-.- b,t '- < '' , ).. b "'- 

b(t 2 -t|) 



\ . i. n -btt r t ) w b(t 2 -t,) -b<t 2 -t,L uM btt 2 -t,). , A1C%1 
X| - k g Ci(l-e v ) (e * ' +e * ' ) ib<t(Ue * . ) (A1B) 
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1-e 

Since the numerator is positive and the denominator is negative, X; ie 
negative. 

Taken together, these lemmas show that acceleration. passes through coast to 
deceleration* There is no possible variation in this scheme. It is also 
possible to show the movement cannot start by coasting followed by 
acceleration. 
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Appendix B 

A surprising limitation on the value of t 2 -tj, the duration of the 
coasting time, has been found. The switching curve during acceleration is: 

X 2 (t)+«+x 2 (t) - X 2 (to)e b(t " to, + ^+k t Cj L (l-e b(t " t ° ) ) + o (Bl) 

b 
. -b(t-t ). 
+ k-ci (1-e v ) 

"V 
The slope of this switching curve is: 

b(t-tn) ,.» ,. » » . » , -b(t-tn) ,__, 

e (bX 2 (t )-Xi-k,C|) + k,c,e ° (B2> 

At the first switching time t,, the switching function (Bl) is zero. 
Rearranging (Bl) for t«tj, 

e b(t,-t ) (X9(tn) . ( Xl+ktC|) ) . . ( x.^c.) - a - k f c 1 (l-e" b(tl " t ° ) ) (B3) 

b b b 

Substituting (B3) into <B2), the slops at t, is: 

- 2k,c,(l-e~ b(t, " to) ) - Xj - bo (B4) 

Substituting for X| from (A16), 

>.c,a-.- b "'-'°' l .2«-'" < '-'- , -. b "'-'' ) ) t 2b. CBS) 

Since slope(tj)>8 and since the denominator is positive, so is the 
numerator. 



. „ -b<t r t ).,„ -b(t 2 -t.) b(t 2 -t,), _ 
k B cj(l-e ' ° )(2+e z ' -e z ' ) + 2ba i 8 



(B6) 
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„ -b(t 2 -t|) b(t 2 -t,) , 
Because 2+e -e is a decreasing function of t 2 , at some 

point (B6) becomes zero. Solving then for e 2 ~ ' j 

b. ♦ VdMtk.cd-e' 1 '" 1 '''' 1 ))' + V.(l-' b " l ' t<, V + 1 (B7» 
K,c,U-.- b,, '- , '> , > 

As t 2 increases, t| Mill decrease. Honever, tj does not decrease enough to 
offset the effect of the t 2 increase. If a«0» (B7) reduces to 

e b(t2 " t » ) - 1 + 4z (B7a) 

Strangely, in this circumstance t 2 -tj depends only on b. 
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